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Abstract 



We present a simple exact analytical solution, using the Weyl-Titchmarsh-Kodaira spectral 
theorem, for the spectral function of the one-dimensional diatomic molecule model consisting 
of two attractive delta function wells in the presence of a static external electric field. For 
sufficiently deep and far apart wells, this molecule supports both an even and an odd state, and 
the introduction of a static electric field turns these bound states into quasi-bound states which 
are Stark shifted and broadened. The continuum spectrum also inherits an intricate pattern of 
resonances which reflect the competition between resonant scattering between the two atomic 
wells and between the linear potential and one or both atomic well(s). All results are analytic 
and can be easily plotted. The relation between the large orders of the divergent perturbative 
Stark shift series and the non-perturbative widths of quasi-bound levels is studied. 



1 Introduction 



The electronic structure of atoms and molecules is usefully probed by external electric and magnetic 
fields. The development of intense lasers has permitted such probing in regimes where simple per- 
turbative treatments are not valid and one must use a nonperturbative semiclassical approximation 
or a numerical approach. Experiments with molecules display a much richer range of phenomena 
than with atoms, due to the additional molecular degrees of freedom These include above 

threshold ionization jl| 1^1 , multiple ionization ^ [7\ , alignment effects [S] , electron localization , 
non-sequential double ionization ^H] , direct excitation , stabilization > dissociative recombi- 
nation J3]) and separation effects ^3]. However, molecules (even the simplest diatomic molecules) 
and their ionization are clearly more difficult to treat theoretically. Many approximation techniques 
have been developed and applied to atomic ionization processes jElElElElEj) but much less 
is known for molecular systems. Realistic calculations are rather complicated and one loses some 
of the physical intuition that can often be gained from simple models. In this paper we present 
the exact analytical solution for a simple molecular ionization model. The molecule is taken to be 
one-dimensional. This approximation is remarkably good in the strong-field regime where the ion- 
ization is predominantly along the field direction, so that the system is effectively one-dimensional 
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[S]. The simplest such one-dimensional molecule consists of two atomic wells^ represented by at- 
tractive delta function wells of strength g, separated by a distance 2a. This is a well-known soluble 
model |2L)j . This molecule always supports an even parity bound ground state and a continuum, 
and if ag > 1 it also supports an odd parity bound excited state. There is a long tradition of using 
model potentials such as zero-range potentials in atomic and molecular physics j21j . The presence 
of an external electrostatic field, of field strength F, dramatically changes the basic character of 
the spectrum, converting the bound states into quasi-bound states and modifying the resonance 
structure of the continuum. These spectral changes are seen directly in the spectral function p{E) : 
the quasi-bound states are poles of p{E) at complex values of the energy E, where the real part of 
the pole gives the energy location of the quasi-bound level and the imaginary part gives the width, 
and hence lifetime, of the level. Since this molecular ionization model is exactly soluble, we can 
easily investigate the dependence of these quasi-bound levels on the relevant physical parameters 
- the field strength F, the atomic well depth g, and the atomic separation parameter o. The 
same applies for the " continuum" , where resonance structures appear due to the delicate interplay 
between tunneling, binding and scattering effects. 

The solution presented here uses the Weyl-Titchmarsh-Kodaira (WTK) spectral theorem j221 1231 
I24| I25j . This spectral theorem expresses the completeness of the wavefunctions of the Schrodinger 
equation in a general way that applies not just to the familiar discrete spectrum models (such as 
the infinite square well or the harmonic oscillator), but also to systems with discrete and continuum 
spectra (such as the finite square well or the hydrogen atom), and even to systems with a purely 
continuous spectrum, such as for ionization problems where there are no true bound states. The 
WTK approach is well suited for numerical implementation and has been applied long ago to 
the Stark effect in atomic hydrogen |26j. An interesting soluble one-dimensional atomic model 
consisting of a single finite square well on the half-line is solved using the WTK method in j2Zl • More 
recently, the numerical WTK approach has been used for studying strong-field ionization effects 
in effectively one-dimensional diatomic molecules |28| I14j . Various numerical and approximate 
methods for computing resonance locations and widths are compared in j^H]- This current paper is 
complementary to |28j . but the choice to represent the atomic wells by delta function wells makes 
the entire molecular ionization problem analytically solvable, thereby bypassing the numerical part 
of the computation. We also mention that this one-dimensional molecular ionization model has 
been studied in [301 using the solution of the associated Lippman-Schwinger equation. 

This paper is organized as follows. Section 2 contains a summary of the implementation of the 
Weyl-Titchmarsh-Kodaira (WTK) method for computing the spectral function. This is applied 
in Section 3 to the one-dimensional molecular model without an applied electric field. In Section 
4 the electric field is applied and the exact solution for the spectral function is derived. The 
dependence of the spectral function on the various physical parameters is explored through plots 
and also analytically. Section 5 is devoted to the simpler case of atomic ionization, obtained from 
the molecular solution by taking the atomic separation parameter, a, to zero. In this case we 
also study the large orders of perturbation theory for the Stark shift and show its relation to the 
nonperturbative level width. The final Section contains some concluding comments. 

^It is straightforward to generalize the exact solution to the case of unequal atomic well depths, but for simplicity 
here we consider the atomic well depths to be equal. 
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2 Weyl-Titchmarsh-Kodaira Method 



The Weyl-Titchmarsh-Kodaira (WTK) spectral theorem j22| \'2'A\ I24| I25j for quantum mechanical 
Hamiltonians is very general, covering not just simple Hamiltonians like the harmonic oscillator 
which have only bound states, but also Hamiltonians with both bound and continuum states. It 
also extends to ionization problems where the spectrum is purely continuum. Indeed, in his classic 
book 123], Titchmarsh solves the half- line "atomic" problem of a binding delta function potential 
well plus a constant electric field. This approach can be applied directly to any one-dimensional or 
radial Schrodinger problem. The WTK method can be summarized as follows. 

Consider the Schrodinger equation (we work in units where ^ = 1)^ 

-ip"{x) + V{x)ip = EiP{x) (1) 

on the real line x G (— oo,+cxd). Pick some point, chosen without loss of generality to be x = 0, 
and normalize the two independent solutions of u{x,E) and v{x,E), so that their Wronskian 
is equal to 1 at that point by choosing: 

u{0,E) = l , u'{0,E)=0 

v{0,E)=0 , v'{0,E) = -l (2) 

Next, integrate (numerically or analytically) each of these solutions out towards x = +oo and 
X = — oo, producing four functions u± (x, E), and v± (x, E). The WTK method involves constructing 
particular linear combinations of these functions such that these combinations are normalizable 
on the intervals (0, +oo) and (— oo,0), when the energy E has a small positive imaginary part: 
E ^ E -\- ie. Specifically, going towards the right, we construct 

i>+{x,E) = u+{x,E) + m+{E)v+{x,E) (3) 

such that it is normalizable on the interval (0, +oo) when the energy E has a small positive imagi- 
nary part. This determines the coefficient function mj^{E). Similarly, going to the left, we construct 
the linear combination 



ii)-{x,E) = u-{x,E)+m^{E)v^{x,E) (4) 

such that it is normalizable on the interval (— oo,0) when the energy E has a small positive imag- 
inary part. This determines the coefficient function m-{E). The spectral function, and the com- 
pleteness of the wavefunctions, can be expressed in terms of these coefficient functions m±[E) 
[22» i22< .24, .2^ . For example, the spectral function is 

piE) = hm '-Im (-^.(E + ^^)^-iEye) + l 
e^OTT \ m+{E + ie) - m^{E + ie) 

This WTK method is well suited to numerical implementation for arbitrary potential wells j28j . 
In the models studied in this paper the situation is even simpler, since the independent solutions 
u{x,E) and v{x,E) are known in analytic form for all x; the "integration" process simply involves 
applying the correct continuity and discontinuity boundary conditions at the locations of the two 
delta function potentials. This will be shown in detail in Section |31 
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3 One-Dimensional Molecular Model 



We first review the model without the electric field. We choose the following simple double-delta- 
function potential to represent the diatomic molecule: 

V{x) = -g [d{x + a) + d{x - a)] (6) 

where g > 0. This potential has two binding delta-function wells located at x = ^a, and chosen for 
simplicity to have equal strength —g. This problem, without an external electric field, is a standard 
problem in quantum mechanics courses [20]. The potential © always supports a bound ground 
state which has even parity, and if ga > 1 it also supports a bound excited state which has odd 
parity. In this section we present the WTK solution as an introductory illustration of the method. 



In the vicinity of a; = the potential © vanishes. Thus, the two independent solutions, u{x, E) 
and v{x,E), of the Schrodinger equation which satisfy the normalization conditions @ are 

u{x,E) = cos{VEx) , v{x,E) = (7) 

v E 



Integrating to the right, these solutions remain valid until we reach the right-hand delta- function 
at X = +a, at which point we impose the standard boundary conditions pUj : 



/ / \ / / N dip 
^(a - e) = -0(0 + e) , — 

ax 



= -g^{a) (8) 

These conditions determine the solutions for x > a to be: 

. , /-= ^ „ sin(v^x) , , 

u+{x,E) = Acos{^x) + B — ' (9) 



E 

v+{x, E)=C cos{VEx) + D ^™(^^) (10) 

V E 

where the coefficients are 

A = l + - ^ J. ^ , B = -gcos^(aVE) (11) 

C = -%sin\aVE) , = -1 + ^ ""^'"^^ (12) 
E 2Ve 

When the energy E has a small positive imaginary part, the linear combination -0+ = -|- rri-^v^ 
in (jSJ will be normalizable on (0, +co) if the exp(— i\/^x) part is eliminated. This determines the 
function m+{E) in ^ to be 

VEA + iB 



Vec + zD 

iE - gy/E cos{ay/E) e^"^ 
'E - g sni{ay/E) e*"^^ 



(13) 
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Figure 1: Plots of the spectral function for the free one-dimensional molecular potential ©• For 
ag < 1 there is only one bound state, but for ag > 1 there are two bound states. 

Since the potential (jEJ is symmetric, it follows that m^{E) = —m-\.{E). Thus, the spectral function 
© is determined to be: 

piE) = \\m.—Im{ -m+{E + ie) ^ -) (14) 

Bound states appear on the negative real energy axis as poles of the spectral function dJ. The 
pole of m_|_, satisfying the transcendental equation 

l + = (15) 

9 

corresponds to the even parity bound state. From 1)141) . the zero of m+, satisfying the transcendental 
equation 




also gives a pole of the spectral function, and corresponds to the odd parity bound state (if it 
exists). There is always an even bound state, given by the pole of m+{E) satisfying (jlSj) . ga > 1 
there is also an odd bound state, given by the zero of m^(E) satisfying (fT6|) . Some plots of the 
spectral function ((TUl are shown in Figs. ^andEJ Note the appearance of the single bound state 
when ga < 1, but of two bound states when ga > 1. In Fig. ^the width of the bound state peaks 
in these plots is artificial, as we have kept the small imaginary part e of the energy nonzero in order 
to show the peaks. In the true e — > limit these bound state peaks have zero width, and so would 
not show up on the plot. Fig. |21 shows the continuum part of the spectrum - notice the periodic 
behavior of the spectral function, due to resonant backscattering between the two delta- functions. 
This periodicity is determined by the separation 2a between the wells, and the well strength g. 
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Figure 2: Plots of the continuum part of the spectral function for the free one-dimensional molecular 
potential Note the periodic behavior due to resonant back-scattering between the two wells. 
The left-hand plot is for ag = 0.5, while the right-hand plot has ag = 1.5. 



4 One Dimensional Molecular Model With an Electric Field 



In this section we solve the problem of the one-dimensional molecular potential © studied in 

the previous section, with an additional external static electric field, of magnitude F. Thus, the 
potential is 

V{x) = -g [6{x + a) + 6{x - a)] - F x (17) 

where g > 0, and F > 0. The potential H17() has no bound states when F ^ 0, but it does have 
quasi-bound states. The corresponding Schrodinger equation ^ is analytically soluble since the 
linearly independent solutions are Airy functions. In fact, the WTK solution proceeds exactly as in 
the case without the electric field, except that the basic trigonometric solutions in ((Jj) are replaced 
by Airy functions. 



4.1 WTK Solution for the Spectral Function 



In the vicinity of x = 0, the independent solutions satisfying the Wronskian normalization condition 
© are 



(Fx + E) 



JP2/3 

where Ai and Bi are Airy functions |31j, and the coefficients needed to satisfy the normalization 
conditions (j^J are 

A^-He) = nB^' (-^) , B^^\e) = -.A^ (-^) 
A^■"\E) = -.F-y^B^{-^^ , B^^\e)=.F^^'^m[-^^ (19) 
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Here we have made use of the fundamental Airy function Wronskian identity j31j : 

Ai{x)Bi'{x) - Ai{x)Bi{x) = - , V x G M (20) 

vr 

As before, integrating to the right, the solutions in (|18|) remain valid until we reach the right-hand 
delta-function well at rc = +a, at which point we apply the delta- function boundary conditions 
This determines the solutions in the region x > a to be: 



p2/3 J ' + " \ p2/3 



where the coefficients are 



Similarly, integrating to the left, the solutions 1)18^ remain valid until we reach the left-hand delta- 
function well at X = —a, at which point we apply the boundary conditions Q. This determines 
the solutions in the region x < —a to be: 

._(..^) - A<,>Ai{J^±^).B<,>Bi[J-^) (23) 

where the coefficients are 

E 



F2/3 



F2/3 
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= g-n'^F-^/^Bi 




-gir^F- 






y ti ± 


-^/^Bi 




+gTT^F- 


-^I^Bi 




= -g^^F' 


'^I^Bi 




+gTT^F~ 


-''I^Ai 



(E-Fa)\^JJE Fa) ,, 



p2/3 



{E-Fa) y , 

p2/3 ) \ F2/3 j 



p2/3 



E \ 



F2/3 J 

E \ 



E 



F2/3 PV3 p2/3) 

^72/3 ) \ F2/3 ) \ F2/3 ) 

p2/3 J y p2/3 J 



p2/3 

{E - Fa) 
F273 



E 

>2/3 



(24) 



Thus, the independent solutions to the Schrodinger equation for the potential H17jl are known 
analytically for all x. As described in Section |2 the WTK method involves finding the functions 
m±{E) such that the linear combinations = u± + m±v± are normalizable as x — > ±00, when 
the energy has a small positive imaginary part. 



In the region x > a, the normalizable solution has the form 

( {Fx + E) \ ( (Fx + E) 

1P+[X, E) (X Al [ ;:;:r7:^ — I — I Bl ' 



This determines m^{E) to be 



p2/3 



B^:;'\e) + iA']^>{E) 



p2/3 



m+{E) 



B^::^\E)+tA^l\E) 

In the region x < —a, the normalizable solution has the form 

{Fx + E)' 



'4>-{x, E) oc Ai 



This determines m-[E) to be 



m_{E) 



p2/3 



b^^\e) 

'WiE) 



(25) 



(26) 



(27) 



(28) 



Given these expressions for m±{E), the spectral function is given by ©, which we can write in 
terms of the coefficients ()22|) and 1)24^ as 



p(E) = lim —Im 

IT 







(E+ie)- 




(E- 


Ne) [b^"' {E+ie)+tA'-_^^ (E+ie) 






iE+ie)A 






-ie) [b^"' {E+ie)+iA''_^'' (E+ie) 



(29) 



This is an analytic expression for the exact spectral function for the potential (|17() . In the remainder 
of this Section we discuss the physical properties of this spectral function, using plots and analytical 
methods. 
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4.2 Plots of the Spectral Function 



Before discussing the analytic properties of the spectral function H29() . we present some plots which 
illustrate how the spectral function depends on the physical parameters, in order to develop some 
intuition for the physical processes involved. In addition, some animations showing how the spectral 
function changes, both in the quasi-bound and in the quasi-continuum region, as we vary the field 
strength F, the atomic separation parameter a, or the atomic well depth parameter g can be found 
at 



4.2.1 Dependence on the electric field strength F. 



p(E) p(E) 




Figure 3: Plots of the spectral function (|29() [solid lines] illustrating how the quasi-bound states 
change as the electric field strength F varies. The dashed line shows the corresponding F = 
molecular spectrum (jl4l In these plots, the atomic well strength is g = 1.5, and the separation 
parameter is a = 1, so the free molecule has two bound states. Notice that as F increases these 
two states are Stark shifted and broadened. 

In the absence of the electric field, the molecule has one or two bound states depending on the 
well separation parameter, a, and well strength, g. And the continuum exhibits structure due to 
Ramsauer-Townsend resonances in the scattering between the two delta wells. These two features 
are illustrated clearly in Figures 1 and 2. When the field strength F is nonzero, the bound states 
become quasi-bound states with a nonzero width, and their central values are Stark shifted. These 
effects are illustrated in Figure 3, which shows a molecule with two bound states subjected to an 



9 




p(E) p(E) 




1?0^ 



Figure 4: Plots of the spectral function ()29() [solid lines] illustrating how the continuum part of the 
spectrum changes as the electric field strength F varies. The dashed line shows the corresponding 
F = molecular spectrum dJ. In these plots, the atomic well strength is g = 1.5, and the 
separation parameter is a = 1. Notice the two oscillation scales - the longer one is set by the free 
case [dashed line], while the rapid oscillation is set by the field strength F, with the free spectral 
function providing an average. As F decreases in magnitude, the oscillations become more rapid, 
eventually averaging out to the free case. 



external electric field of strengths F = 0.05,0.1,0.5,1. Note that the quasi-bound states broaden 
as the field increases, with the higher state being broader since it is less deeply bound and so can 
tunnel more easily. The quasi-bound states are Stark shifted in opposite directions: the lower state 
is Stark shifted down in energy, while the higher state is Stark shifted up in energy. The effect 
of the external electrostatic field on the continuum states is shown in Figure 4. In these plots 
the dashed line shows the continuum spectral function when F = 0. The solid lines show the 
spectral function for various values of F: F = 5, 10, 15, 20. Note that the free spectral function 
provides an average for the F > spectral function. The scale of the oscillations of this average 
function is clearly independent of F, being determined by the combination ag, as illustrated in 
Figure 2. For F > there is an additional scale in the "continuum" spectrum, and as F decreases 
to zero the rapid oscillations become more and more rapid, and eventually average out to the free 
spectral function. As F increases the period of this oscillation increases. These oscillations can be 
correlated approximately with the zeros of the Airy function Ai{—E/F'^^^), since these zeros give 
the energies of the half-wedge potential well which has V = —Fx for x < 0, with an infinite barrier 
at X = 0. This is illustrated in Figures 5 and 6, for a = and a = 1, respectively. Note that the 
peaks of the spectral function (the solid line) coincide roughly with the zeros of Ai{-E/F^/^), the 



10 




Figure 5: Plots of the spectral function ()29() [solid lines] illustrating how the continuum part of the 
spectrum correlates with the zeros of the Airy function Ai{—E/F'^^^) [dashed line], which are the 
energies of the wedge potential V = —Fx with an infinite wall at some point. This is an illustration 
of the Ramsauer effect. These plots are for a = and g = 1.5, and various values of the electric 
field strength F, as shown. 

dashed line. The agreement is quite good, even for the molecular model having a = 1. This is an 
example of Ramsauer- Townsend resonance, with the electron backscattering off the two delta wells 
providing the large period oscillations, and the electron scattering off the delta wells and the linear 
electrostatic potential providing the shorter period oscillations. 



4.2.2 Dependence on the atomic well separation parameter, a. 

The dependence of the spectral function (|29j) on the atomic well separation parameter, o, is illus- 
trated in Figures 7 and 8, for the "bound" and "continuum" parts of the spectrum, respectively. 
In the zero field case the spectral function has two bound states ii ga> 1. These are shown as the 
dashed curves in Figure 7. As the separation increases these two bound states approach the same 
energy, becoming degenerate in the limit a ^ oo. This is because in the large separation limit the 
tunneling which mixes the two levels becomes suppressed and the two atoms become essentially 
independent of one another. Thus the bound state spectrum approaches that of a single atomic 
well. On the other hand, if the field strength is nonzero the quasi-bound states do not become 
degenerate in the limit of large separation. This is shown by the solid curves in Figure 7. Instead, 
the two states move away from each other; the even state is Stark shifted further down in energy 
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Figure 6: Plots of the spectral function ()29() [solid lines] illustrating how the continuum part of the 
spectrum correlates with the zeros of the Airy function Ai{—E/F'^^^) [dashed line], which are the 
energies of the wedge potential V = —Fx with an infinite wall at some point. This is an illustration 
of the Ramsauer effect. These plots are for a = 1 and g = 1.5, and various values of the electric 
field strength F, as shown. 

and the odd state is Stark shifted further up in energy. This is because for large a the tunneling is 
essentially from each well independently, with one Stark shifted up and the other down, depending 
on the parity of the original state. These quasi-bound levels also become narrower in width as 
there is a larger barrier through which the electron must tunnel as a increases. The effect on the 
"continuum" part of the spectrum is shown in Figure 8. Here F = 0.15 and the spectrum exhibits 
resonances which oscillate about the -F = case, which is shown by the dashed lines. As the 
well separation changes, the average function changes, due to the Ramsauer- Townsend resonance 
between the two wells. Note that for nonzero field strength the spectral function oscillates rapidly 
with energy, but still follows the free field average. Thus, the F > spectral function represents a 
competition between the resonance between the two wells, which depends strongly on a, and the 
resonance between a given well and the linear potential, which is not sensitive to a. 



4.2.3 Dependence on the atomic well strength, g. 



The dependence of the spectral function (|29j) on the atomic well depth, g, is illustrated in Figures 
9 and 10, for the "bound" and "continuum" parts of the spectrum, respectively. In the -F = case, 
the difference in energy between the even and odd states gets smaller as the well strength increases, 
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Figure 7: Plots of the spectral function [solid lines] illustrating how the quasi-bound part 
of the spectrum changes as the well separation parameter a varies. The dashed line shows the 
corresponding F = molecular spectrum ()14|) . For the solid lines, the field strength is -F = 0.1, 
the well strength is = 1.5, and a ranges through 1.5, 2, 2.5 and 3. Notice that as a increases the 
quasi-bound states become narrower in width and they move apart. In contrast, the free bound 
states move together, eventually becoming degenerate as a ^ oo. 

which is similar to what happens when the well separation increases. This is because the tunneling 
mixing is more highly suppressed as the states become more deeply bound. The difference between 
the limit of large well separation and the limit of large well strength is that as the well separation 
increases the two bounded states approach the same energy, whereas when the well strength is 
increased the states become degenerate but their energies tend to — oo as ^ oo. However, if the 
electric field is applied, the two states do not become degenerate in the large well strength limit but 
instead keep a nonzero relative distance between each other. The distance between the two states 
approaches 2Fa as well strength increases, which is easily explained by the following argument. 
If the well strength is very strong we can think of a particle being localized at a single well. If a 
uniform electric field of strength F pointing in the positive x-direction is applied, a particle localized 
at the left well will increase in energy by Fa while the energy of a particle localized at the right 
well will decrease by Fa. Since the unperturbed states have the same energy the energy difference 
is just 2Fa. This can be seen in the last panels of Figure 9, for which F = 0.1 and a = 1.5, so 
2Fa = 0.3, which is roughly the separation between the two quasi-bound levels. Also note that the 
quasi-bound state peaks become narrower as g increases, as the levels are more deeply bound. 

Backscattering resonances in the continuum become more prominent as the well strength increases. 
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Figure 8: Plots of the spectral function (|29() [solid lines] illustrating how the continuum part 
of the spectrum changes as the well separation parameter a varies. The dashed line shows the 
corresponding F = molecular spectrum ()14() . In these plots, the field strength is F = 0.15, 
and the well strength is g = 1. Notice that as a increases a second quasi-bound state peels off the 
positive energy continuum and forms a Stark-shifted pair around the two F = bound states which 
exist for ag > 1. Also note that in the F = case, as a increases the splitting between the two 
bound states becomes vanishingly small, as the tunneling between the two wells is suppressed. In 
the continuum, the average function varies with a due to resonances in the backscattering between 
the two wells. Even with nonzero F the spectral function follows this average closely as a varies. 

because the resonances are sharper, since the scattering potentials are deeper. If the field is weak, 
the WKB approximation for a infinite well potential in a electric field can be used to find the 
energies of the backscattering states in the strong well strength limit. But these peaks also follow 
the free field average (the dashed lines in Figure 10), which is due to the resonance between the 
wells without the electric field. Thus, the F > spectral function represents a competition between 
the resonance between the two wells, which depends strongly on g, and the resonance between a 
given well and the linear potential, which is also sensitive to the well strength g. 



4.3 Analytic Properties of Spectral Function: Stark Shifts and Level Widths. 

With the electric field present there are no true bound states, but there are quasi-bound states. 
The location of these quasi-bound states are given by the real parts of the poles of the spectral 
function (j29]) . and their widths are given by the imaginary part of these poles. These poles are 
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Figure 9: Plots of the spectral function (|29() [solid lines] illustrating how the quasi-bound part of the 
spectrum changes as the well depth parameter g varies. The dashed line shows the corresponding 
F = molecular spectrum ()14() . In these plots, the field strength is -F = 0.1, and the well separation 
parameter is a = 1.5. 
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Figure 10: Plots of the spectral function (|29|) [solid lines] illustrating how the continuum part of the 
spectrum changes as the well depth parameter g varies. The dashed line shows the corresponding 
F = molecular spectrum p4|) . In these plots, the field strength is -F = 0.5, and the well separation 
parameter is a = 1. Notice that as g increases a second quasi-bound state peels off the positive 
energy continuum and forms a Stark-shifted pair around the two F = bound states which exist 
for ag > 1. 
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given by the zeros of: 



(E + Fa) 



X Ai 



Bi 



{E - Fa) 

^273 

{E - Fa) 



Ai 



{E + Fa) 

p2/3 



iBi 



(E + FaY 
p2/3 

{E + Fa) 



F2/3 



-girF' 



-2/3 



^2/3 / 



+ Ai 



{E - Fa) 

F273 



Bi 



{E-Fa) \ 
F^ ) 



p2/3 



p2/3 



{E + Fa) 



(30) 



In the weak field hmit we can use the fohowing asymptotic expansions [HT] of the Airy functions to 
find perturbative solutions for the zeros of H30|): 



Ai{z) 
Bi{z) 



A:=0 

oo 



1/4 X] 



fc=0 ^ 



Cfc 



(31) 
(32) 



where C, = 1^^''^, and the expansion coefficients Ck are 



Cfe 



r(3A; + 1/2) 
W'k\T{k + 1/2) 



(33) 



Clearly, only even powers of F will appear in the perturbative expansion for the real parts of the 
quasi-energies. So, we define the expansion 



E 



oo 
n=0 



2n 



(34) 



To find the real parts of the quasi-energies, we can ignore the imaginary parts of (|3()j) . which are 
anyway exponentially suppressed in the weak field limit. It is a straightforward exercise to expand 
(|Hn|) in powers of the field strength F, using for example Mathematica jSS]- The leading order F^ 
term produces the equation 



4an 



4^+ (l-e-^^sV^ 







(35) 



whose solutions are just the solutions of the transcendental equations ((T3|) and l(TB|) derived in the 
previous section for the bound states in the free field case. A solution to satisfying (|T^ is 
an even bound state of the F = potential while a solution to satisfying ((TB|) is an odd 
bound state of ©• 

The first correction to these bound states comes from the F'^ term in the expansion of H3U|) . which 
leads to the following expression for ai in terms of cq: 

^ _ 15(1 - 2^^) - 15ag(l - 2V^)2 + 12(ag)2ao(l - 4^) - A{ag)^aoil - 2^)^ 

48a2(l-2V^)(l-a<7(l-2^)) ^ ^ 
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Thus, to first nontrivial order, the Stark shifted energy is 

E = -aog'^ - + . . . (37) 

To find the shift for the even bound state, we find the solution oq of (|52j) which also satisfies (|15|). 
and then insert this value of oq into H36() to find the corresponding ai. It is straightforward to 
continue this to higher orders. (We study details of the higher orders in the Section 1^1 for the 
special "atomic" case where the atomic separation parameter a = 0.) For the odd bound state, we 
must first find if there is such a solution for uq to (|52|) satisfying (|16|) . This odd solution will exist 
if ga > 1. If it exists, then the corresponding Stark shift is obtained by inserting this value of oq 
into the expression (j36|) for ai. 

The widths of the quasi-bound states can be derived from the imaginary part of the poles of the 
spectral function. Because of the complicated dependence of the quasi-bound state energies on the 
system parameters F, g and a, it is difficult to derive simple analytical expressions for the line 
widths. However, in the large atomic separation limit, as a — > oo, we can use WKB to approximate 
the tunneling rate, and hence the line width, as 



r ~ g\/ Eq — Fa exp 



i{Eo-Faf/^' 



(38) 



3F 

where E = —Eq — iV, and —Eq is the full Stark-shifted energy of the lower quasi-bound state. 

In the limit of infinitely large well strength, g ^ oo, the g^ term in ()3U() dominates and the zeros 
of ()3U() lie on the positive real axis. If the field strength F is small, the first order asymptotic 
expansion of the Airy functions |31j : 

M-^) ^sm(c+^) (39) 



^zV4 4 



for Re{z) ^ can be used to estimate the location of the zeros of (|30p . These approximate zeros 
are determined by the expressions: 

2(^-Fa)3/2 / 1\ 

- '"'--Itt (41) 



3 F V 4 

2 (E + Faf/^ - (E - Faf/^ , , 
'- ^ '- = nvr (42) 

3 F ^ ' 

Note that expression (|4H) approximates the zeros of Ai ^— ^^2/3"^ ^ while (|42|) approximates the 

zeros of Ai ^— ^^^"^ ^ Bi ^— ''^2/3°'* ^ ~ ^— ^^2/3"^ ^ ^— ^^^"^ ^ Solutions for expression 
H41|) are the energies obtained using the WKB approximation for a half wedge potential with an 
infinitely high wall at x = —a. Therefore, we interpret the energies satisfying H41|) as those of 
the backscattering states between the left well and the electric field, which is the Ramsauer effect 
illustrated in Figures 5 and 6. Similarly, expression (|42j) yields the same set of energies obtained 
from the WKB approximation for a potential well in an electric field with infinitely high walls at 
X = —a and x = a, and therefore (|42|) determines the energies of backscattering states between 
the two wells, also in the presence of the electric field. In general, the exact resonances reflect a 
competition of these scatterings amongst the linear potential and the atomic wells. 
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5 Single Well Potential With an Electric Field 



A nice feature of our molecular analysis is that the "atomic" analogue of the molecular model 
studied in Section |1] can be obtained simply by setting the separation parameter, a, of the two wells 
to zero. All the expressions carry over smoothly in this a — > limit. The corresponding potential 
is (note that the delta function strength becomes 2g in this limit): 

V{x) = -2g6{x) - Fx (43) 

This atomic problem has also been discussed in terms of the corresponding Green's functions in 

[HSIISI. 



5.1 WTK Solution for the Atomic Spectral Function 



The spectral function is given, as before in (|29j) . by 



p{E) = lim —Im 

e— >0 TT 



(E+ie) ^B'_;;:'{E+ie)+iA'_^' (E+ie)^ -B^"' {E+ie) \ B^"-' {E+ie)+iA'^' (E+ie) 

where the coefficient functions in H22|) and (|24() now simplify to : 
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Taylor expanding in e leads to a more explicit expression for the spectral function 

1 
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(44) 



(45) 
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x{2gF-'^/^Ai 



F2/3 J > 



E 

>2/3 



E 

>2/3 



— Y 



(46) 



5.2 Large-Order Perturbation Theory and Level Widths. 



It is well known that perturbation theory for an unstable quasi-bound state produces a divergent 
non-alternating series, with which an imaginary part can be associated using Borel techniques. 
This was first explicitly investigated for the unstable + gx^ anharmonic potential ;36j , then in 
extensive detail for the x'^ + Xx^ potential 37^, and is a very general property of perturbation theory 
[SHI EH]- The same ideas apply to the Stark effect problem, as has been investigated exhaustively 
for atomic systems (lUl HTl EH] . 



The quasi-bound states correspond to poles of the spectral function, namely solutions to 

2 



(- " VI 


2 


V FV-i) 


+ 



2TTgF-^/^Bi 



E 

>2/3 



Ai 



E 

>2/3 ; 







(47) 



The real part of the quasi-energy can be found by making a perturbative expansion for the real 
part, Sreab of the energy 



2n 



(48) 



as in (|.S4|) . When the field strength F vanishes, there is a single (even) bound state at E = —g^. 
Thus, ao = 1, as is consistent with when a — > 0. The width of the quasi-bound state can be 
estimated by writing 



(49) 



where we expect -Eimag to be exponentially small. Indeed, expanding the imaginary part of (|47|) 
immediately leads to the leading behavior 



E- ~ 



-g^ exp 



3F 



(50) 



in agreement with the a — > limit of the molecular case H38() . 



We now show how this is consistent with an analysis of the divergence of the perturbative expan- 
sion for the real part of the quasi-bound energy level. First, note that the Ai"^ term in H47|) is 
exponentially small in the small F limit, and so can be neglected. Thus the real part is determined 
by 
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n 



NlaJ 



N 



J (lead) 



R 



(lead) 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



1 

_5_ 

1 

64 
10625 



Iff 



366fW625 



3276 
9378457 



!125 



131072 
269028257953125 
4194304 
129011616275390625 

16777216 
159621687625662109375 
134217728 
123839968932138228515625 



471147487' 



j?i?^?246943359375 



5392128838§t^§^iM^10673828125 

34359738368 
1462185114846262625626556396484375 

1448716002754jl§J^??f?991 76025390625 



679691840600372; 



mmm 



i225469970703125 



70368744177664 
142608185435906164633493702703533111572265625 
281474976710656 
679265718819054465192747030828993319061279296875 



18224958526834818420 



2« 



§35^639728546142578125 



218881880317532293574 



mi 



3254740992 

595827650045023616790771484375 



144115188075855872 
73105891881984796538857909985635411709301443950653076171875 



1 

0.3125 
0.859 
10.376 
263.214 
11198.1 
715520. 
6.4.10'^ 
7.7. 10^ 



576460752303423488 



1.2 

2.3 
5.5 
1.6 
5.3 
2.1 
9.7 
5.1 
3.0 
2.0 
1.5 
1.3 



1012 

10^4 
10^6 

1019 

1021 
1024 
1026 

1029 
10^2 
1035 
1038 

10 
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0.872665 
0.71106 
0.763133 
0.819423 
0.860776 
0.888895 
0.908071 
0.921614 
0.931602 
0.939271 
0.945356 
0.950311 
0.95443 
0.957911 
0.960894 
0.963479 
0.965742 
0.96774 
0.969516 
0.971107 



1.04760 
0.95104 
0.9533 
0.98733 
1.006199 
1.008710 
1.005397 
1.002243 
1.000556 
0.99992 
0.99978 
0.99981 
0.99987 
0.99991 
0.99995 
0.99997 



Table 1: This table lists the coefficients o™ appearing in the perturbative expansion (|48|) for the 
quasi-energy level fo the atomic system in an electric field. Note that these coefficients are non- 
alternating in sign, and that their magnitude, A^[a„], grows very fast with the perturbative order 
n, as is shown in the third column. The fourth column shows the ratio of the a„ to the leading 
factorial growth rate in (|52|) . and the fifth column gives the fourth-order Richardson extrapolation 
[H| of this ratio, showing its rapid approach to unity. 
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In the small F limit we can use the asymptotic expansions (|31I32|) of the Airy functions to make 
an expansion of (|51|) in powers of (F/g^)'^, thereby successively determining the coefficients a„. 
It is trivial to program this expansion in Mathematica, and the results for the first 21 expansion 
coefficients are shown in Table 1 |35j . 

Notice that all the an have the same sign, and their magnitude grows very rapidly. In fact, the 
leading growth rate at large orders n of the perturbation series is factorially fast: 

2 / S \ ^" 

^(W)^ r(2n) (52) 



The last two columns of Table 1 show how quickly this leading growth rate is achieved. The final 
column uses fourth order Richardson extrapolation for the ratio of the exact coefficients to the 
leading behavior in (|52|) . and the rapid onset to the leading growth shows that after 10 terms one 
is already well and truly into the asymptotic regime described by (|52j) . 

Since the perturbative series H48|) for the energy is divergent and non-alternating, it is not Borel 
summable [llj. However, one can associate an imaginary part with such a divergent series using 
standard Borel dispersion relations [SHI EZl OH IHSj' This technique provides a bridge between 
perturbative and nonperturbative physics, and is well studied in a wide variety of physics contexts 
[38j . The argument can be motivated and illustrated roughly as follows. Consider the following 
alternating divergent series obtained by an asymptotic expansion of the integral 

TrT^~E(-l)"^'"r(2n) (53) 

n=0 

where c > 0. The integral on the left-hand-side of (|53jl is the Borel sum of the alternating divergent 
series of the right-hand-side. This argument does not work for the corresponding nonalternating 
series, Yl'^=o^^^i'^''^)y ^ would encounter a pole on the t integration axis. The behavior at 
such a pole must be defined by some sort of prescription. The standard principal parts prescription 
associates with this non- alternating divergent series a non-perturbative imaginary part |38l I39j 



\n=0 



\ vr 


1" 


1 ~ - exp 




c 



(54) 



Thus, using the leading growth rate (|52|) for the perturbative coefficients we see that this Borel 
approach associates an imaginary part for the quasi-bound energy level equal to 



Im{E) ~ -g^ exp 



(55) 



This agrees precisely with the imaginary part (|5()|) found directly from the properties of the Airy 
functions, as well as with a simple WKB estimate. 



6 Conclusions 



In this paper we have presented the exact analytic solution for the spectral function for the simple 
one-dimensional molecular ionization model of a diatomic molecule represented by two attractive 
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delta function wells in an external static electric field. The Weyl-Titchmarsh-Kodaira spectral 
theorem provides a simple construction for the spectral function in terms of suitably normalized 
solutions to the Schrodinger equation. In this case these solutions are Airy functions, and the 
spectral function can be expressed in closed form in terms of Airy functions. Thus, the spectral 
function can easily be plotted using a program such as Mathematica 35 . The dependence of the 
spectral function on the relevant physical parameters, the field strength F, the well strength g, 
and the well separation parameter a, is illustrated in Section 4 by a collection of plots. This helps 
develop a body of intuition for the behavior of the quasi-bound states as they are Stark shifted and 
broadened, and also for the resonance structures in the "continuum", which reflect a competition 
between Ramsauer-Townsend resonant scattering between the two atomic wells and between one 
or both atomic well(s) and the linear field potential. The most important extension of this model 
would be to consider the effect of time-dependence in the background electric field, which introduces 
yet another physical scale into the problem El BEl QHl 1^01 • 
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